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Abstract 

Spin dynamics of nanomolecules and nanoclusters are analyzed. The nanosizes of 
these objects make it possible to consider them as single-domain magnets with a large 
total spin, where the motion of the spins of all atoms, composing a nanocluster, occurs 
in a coherent way. Another meaning of coherence in spin dynamics is the coherent spin 
motion of several nanomolecules or nanoclusters. Different approaches for treating 
spin dynamics are compared and the main mechanisms influencing the spin motion 
are studied. Spin dynamics of separate magnetic nanomolecules and nanoclusters are 
investigated, as well as the spin dynamics of the ensembles of these nano-objects. 

1 Introduction 

Magnetic nanomolecules and nanoclusters enjoy many similar properties because of which 
the dynamics of their magnetization can be described by the same type of equations. This 
is why, we consider both these nano-objects together. Of course, there is difference in their 
structure and parameters which we shall take into account and characterize them by the 
appropriate models. The detailed description of general physical properties and applications 
of different magnetic nanoparticles can be found in review articles [1-9]. Here we briefly 
mention those of the properties and parameters that will be necessary for the following 
consideration. 

It is worth stressing that there exist two types of magnetic nanoparticles. One large 
class consists of nanoclusters and nanomolecules, whose magnetic moments are formed by 
electron spins. Another type includes nanomolecules that possess magnetic moments solely 
due to polarized proton spins. Examples are propanediol C 3 H 8 2 , butanol C 4 H 9 OH, and 
ammonia NH 3 . In such nanomolecules, there is no any other magnetic moment except that 
caused by polarized protons. So, here the proton magnetic moment is not a contribution, 
but the main object. 

The magnetic moment of an atom is composed of electron and proton moments, with the 
electron magnetic moment \i e = —g^sS = h-y e S = —\ib and the proton magnetic moment 



— 9pI^nS = h'fpS, where g e = 2 and g p = 5.586 are the electron and proton Lande 
factors, fx b = \e\h/2m e and hn = \e\H/2m p are the Bohr and nuclear magnetons, 7 e and 7 P 
are the electron and proton gyromagnetic ratios. Since the proton mass is larger than that 
of an electron, m p /m e ~ 10 3 , the proton magnetic moment is essentially smaller, fx p / fx e ~ 
10~ 3 . The electron and proton radii are r e ~ 10~ 15 cm and r p ~ lCT 13 cm, respectively. An 
atom is called magnetic, when its total magnetic moment is nonzero. The total spin of a 
magnetic atom can be between 1/2 and S ~ 10, hence, its magnetic moment can be of order 
1/Xb — lOfiB- Atom radii are of order ta ~ 10~ 9 — 10" 8 cm. Examples of magnetic atoms are 
Fe (Iron), Co (Cobalt), Ni (Nickel), Gd (Gadolinium), and Cr (Chromium). 

Magnetic nanomolecules are composed of many magnetic atoms an, as is clear from 
their name, are of the nanometer size. An important property of a magnetic nanomolecule 
is that its total magnetic moment can be treated as being due to an effective total spin. 
Generally, the molecule spin can be directed either up or down, with an energy barrier 
between these directions of order Ea ~ 10 — 100 K. At high temperatures, above the blocking 
temperature Tb ~ 1 — 10 K, a magnetic molecule behaves as a superparamagnetic particle, 
whose spin randomly oscillates between the up and down positions. While below the blocking 
temperature the spin is frozen in one of the directions. 

Magnetic nanoclusters are also made of magnetic atoms that are assembled together in 
a random way. This distinguishes them from magnetic molecules, where atoms are strictly 
connected by chemical bonds. The sizes of nanoclusters can be in the range between 1 nm 
and 100 nm, containing about 100— 10 5 atoms. These values define the coherence radius R co h, 
below which a nanocluster is in a single-domain state and can be treated as a large particle 
with an effective spin. A cluster, with a size larger than R co h, separates into domains with 
opposite magnetizations. Similarly to magnetic molecules at low temperature, the magnetic 
moment of a nanocluster, below the blocking temperature Tb ~ 10 — 100 K, is frozen in one 
of two possible directions. The effective spin of a nanocluster is formed by electron spins 
and can be as large as S ~ 100 — 10 5 . 

The often considered nanoclusters are made of the magnetic atoms of Fe, Ni, and Co. 
They can be made of oxides, such as NiO, Fe 2 3 , NiFe 2 4 or alloys, such as Nd 2 Fe 14 B, 
Pr 2 Fe 14 B, Tb 2 Fe 14 B, DyFe 14 B, Pr 2 Co 14 B, SmiFenTii, S mi Fe 10 V 2 , Sm 2 Fe 17 N 23 , Sm 2 Fei 7 C 2 2, 
Sm 2 Coi7, Sm 2 Co5. To protect nanoclusters from oxidation, one coat them with graphene or 
noble metals, forming the double-component nanoclusters, such as Fe-Au, Co-Au , Co-Ag, 
Co-Cu, Co-Pt, Co-Pd, Ni-Au, Ni-Ag, Ni-Pd, and Mn~-Au. The coating is done be means of 
chemical reactions or laser ablation techniques. The nanoclusters are produced by employing 
thermal decomposition, microemulsion reactions, and thermal spraying. 

Magnetic nanoclusters and nanomolecules find numerous applications, among which we 
can mention magnetic chemistry, biomedical imaging, medical treatment, genetic engineer- 
ing, waste cleaning, information storage, quantum computing, and creation of radiation 
devices. Since both nanomolecules and nanoclusters possess many common properties and 
can be considered as single particles with a large spin, we shall often talk on nanoclusters, 
implying that similar effects can be realized with both of them, molecules as well as clusters. 

The use of these nano-objects requires the existence of two properties that contradict 
each other. From one side, to be able to keep memory, a cluster has to enjoy a stable state 
with its spin frozen in one direction. But from another side, in order to be able to manipulate 
the cluster magnetization, there should exist a way of suppressing the anisotropy. And it 
is necessary that the spin manipulation could be done sufficiently fast, so that the cluster 
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magnetization could be quickly reversed. Recall that thermal reversal is characterized by 
the Arrhenius law giving the longitudinal relaxation time 7\ ~ expl^/Zc^T}, where 
is the anisotropy energy, so that, at temperatures below the blocking temperature, the 
magnetization is frozen. 

Magnetization reversal can be realized by different methods, by applying transverse con- 
stant or alternating magnetic fields and short magnetic field pulses [10]. To achieve fast 
reversal, one needs to find optimal values for the amplitude, frequency, and duration of such 
field pulses. 

A very efficient method of achieving ultrafast magnetization reversal of magnetic nan- 
oclusters has been suggested [11] by employing the acceleration effect caused by a resonator 
feedback field. The efficiency of this method is due to self-optimization of the spin motion 
producing the resonator field acting back on the spins. Historically, this effect was described 
by Purcell [12] and considered by Bloembergen and Pound [13] using classical phenomeno- 
logical equations. Such equations are not sufficient for describing different regimes of spin 
motion. Microscopic theory of spin dynamics has been developed being applied to polarized 
proton spins of such molecules as propanediol C3H8O2, butanol C4H9OH, and ammonia NH 3 
(see review articles [4,14]) and to magnetic molecules [15-19]. 

The aim of the present paper is threefold. First, we concentrate on the spin dynamics of 
nanoclusters, comparing the peculiarity of their spin motion with that of proton and molec- 
ular spins. Second, we analyze the role of other effects, such as the Nyquist-noise triggering 
and Dicke correlation, studying their influence on the spin dynamics of nanoclusters. We 
show that these effects are negligible as compared to the Purcell effect. And, third, we 
compare different approaches to describing spin dynamics, demonstrating the advantage of 
using a microscopic approach based on quantum equations of motion. 



2 Phenomenological classical equations 

Dynamics of the magnetic moment M of a magnetic particle is usually described by the 
classical equation 

^ = -| 7s |MxH e// + R, (1) 

in which 75 is the giromagnetic ratio of the particle with spin S and R is a relaxation term. 
The effective magnetic field is given by the variational derivative H e ff — — 5E/5M of the 
particle energy E. The length of the magnetic moment is conserved, when the right-hand 
side of the equation <iM 2 /dt = 2M • R is zero. 
Choosing the relaxation term in the form 

R=-^Mx(MxH e// ) , (2) 

one gets the Landau-Lifshitz equation, where a is a dissipation parameter and M = |M|. 
Under form (2), |M| is conserved. The equation was initially derived [20] for describing 
energy dissipation in the process of magnetic domain wall motion inside bulk ferromagnetic 
matter. Though it is often applied for treating the dynamics of ferromagnetic particles [21]. 
Taking the relaxation term as 

_ a , , rfM . , 

R = — M x — — , 3 

M dt ' v ' 
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one comes to the Gilbert equation [22]. This equation, up to a renotation of parameters, is 
equivalent to the Landau-Lifshitz equation. Hence, it has the same region of applicability, 
though it is also used for describing the magnetization rotation of magnetic particles [10]. 
Another form of the relaxation term has been advanced by Bloch [23] as 

M x -M* My — M* M z -M; 

J-2 J-2 J l 

where M* is an equilibrium magnetization, e a are unit coordinate vectors, and the relaxation 
parameters are characterized by the longitudinal relaxation time T\ and transverse relaxation 
time T 2 . The latter is also called the dephasing time. For an ensemble of N magnetic particles 
with a large average spin polarization 

1 N 

-sivE^- ( 5 ) 

i=i 

the transverse term has to be renormalized [16,24] as I/T2 = (1 — s 2 ) /T 2 . 

The Landau-Lifshits equation has a single dissipation parameter a and preserves spherical 
symmetry, thus, describing isotropic magnetization rotation. Because of these properties, it 
is appropriate for bulk macroscopic ferromagnetic matter with spherical magnetic symmetry. 
It may also be used for magnetic clusters, possessing this symmetry, which, however, is a 
rather rare case. 

The Bloch equation has two relaxation parameters, T\ and T 2 . Therefore it can describe 
more general situation of anisotropic relaxation, which is more realistic for treating nanoclus- 
ters in a medium or below the blocking temperature, when T 2 < Ti. The Bloch equations 
have been employed for considering the electron and nuclear spin motion in a strongly coher- 
ent regime [25-29] and for spin-polarized 129 Xe gas [30] . But these equations cannot describe 
the whole process of spin relaxation starting from an incoherent quantum stage, for which 
a microscopic approach is necessary [29,31,32]. The initial stage of spin relaxation is trig- 
gered by quantum spin fluctuations that can be identified with nonequilibrium spin waves 
[4,16,31-33]. 



3 Microscopic quantum approach 

In a self-consistent quantum approach, we start with a microscopic spin Hamiltonian H that 
is a functional of spin operators S. The evolution equations are given by the Heisenberg 
equations of motion 

zh |! = [S, H] . (6) 

The advantage of using the quantum approach is in the following. First, it takes into 
account quantum effects that can be important for small clusters. Hence, it is more general. 
Second, at the initial stage of free spin relaxation, quantum spin fluctuations are of principal 
importance, being the triggering mechanism for starting the spin motion. Third, being based 
on an explicit spin Hamiltonian makes it possible to control the used approximations and to 
have well defined system parameters. 
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We assume that a magnetic cluster is inserted into an magnetic coil, of n turns and length 
I, of a resonant electric circuit characterized by resistance R, inductance L, and capacity C . 
The coil axis is taken along the axis x. Moving magnetic moments induce in the coil the 
electric current j described by the Kirchhoff equation 

in which the magnetic flux $ = 4nnM x /cl is formed by the mean transverse magnetization 
M x = fJ-oY2f = i {Sp, where /io = ^7s- Here Ef is an additional electromotive force, if any. 
The resonator natural frequency and circuit damping, respectively, are 

u> = 7 = — . (8) 



The coil current creates the magnetic field 

= ' cl 



B=^j (9) 



dH 
~dt 



that is the solution to the equation 

+ 2 1 H + co 2 f H{t') dt' = ^ , (10) 

Jo dt 

where rj = V/V coi i is the filling factor and 

3=1 

is the transverse magnetization density. The external electromotive force is omitted. The 
field H is the feedback field, created by moving spins and acting back on them. 

4 Dynamics of a single nanocluster 

The typical Hamiltonian of a nanocluster is 

H = -/i B • S - D(S Z ) 2 + D 2 (S X ) 2 + D A [{S x f{S y f + (^(S 2 ) 2 + (S Z ) 2 (S X ) 2 ] , (11) 
where the total magnetic field 

B = B e z + B iex + He x (12) 

consists of an external constant field B , weak transverse anisotropy field B 1 , and the feed- 
back resonator field H. The anisotropy parameters D, D 2 , D A are defined by the particular 
type of considered nanoclusters. 

The main attention will be payed to the investigation of spin dynamics starting from a 
strongly nonequilibrium initial state, where the magnetization is directed opposite to the 
constant external magnetic field B . 
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First, we study the influence of the thermal Nyquist noise of the coil in order to un- 
derstand whether it can trigger the spin motion in a nanocluster. For the thermal-noise 
relaxation time, we find 



tT = — — n — tanh 



2lot 



(13) 



where ut = ksT/h is the thermal frequency defined by temperature T. At low temperatures, 
below the blocking temperature, say at T = 1 K, we have ut ~ 10 12 s _1 . Then the thermal- 
noise relaxation time is 



coil 



UJ 



< 1 



(14) 



On the other side, for the reversal time, caused by the resonator feedback field, we have 

Vcoil 



t ~ 

"rev — 



(15) 



The ratio of the latter to the thermal time (13) is t rev /tT ~ ut/^^S. For the typical values 
T = 1 K, 7 ~ 10 10 s -1 , and S ~ 10 3 , this ratio is small: t rev /tT ~ 10~ 2 . Therefore the 
thermal Nyquist noise does not play any role in the spin dynamics of a nanocluster. 

We have accomplished numerical solution of the evolution equations for nanocluster 
parameters typical of Fe, Ni, and Co nanoclusters. The Zeeman frequency is taken as 
u = 2fi B B /h ~ 10 11 s~ l . For the feedback rate, we have 7 = irqfvjgS /V coi i ~ 10 10 
s _1 . The typical anisotropy parameters satisfy the relations D/(fry ) ~ 10~ 3 , D 2 /(^7o) ~ 
10~ 3 , D 4 /(h-f ) ~ 10~ 10 . At the initial time, the spin is assumed to be directed along the axis 
z. The resonator natural frequency is taken to be in resonance with the Zeeman frequency 
defined by the field B . The behavior of the spin polarization (5) is shown in Fig. 1, where 
we compare the spin motion in the presence of the resonator (h ^ 0) and in the absence of 
the latter (h = 0). Clearly, without the resonator feedback field, the spin is blocked, while 



in the presence of the resonator, it reverses in short time t r 
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Figure 1: Spin reversal of a single nanocluster, 
with parameters typical of nanoclusters made of 
Fe, Ni, and Co. 



5 Dynamics of nanocluster assemblies 

The ensemble of nanoclusters is described by the Hamiltonian 



(16) 



where the indices i,j = 1,2, ... ,N enumerate nanoclusters. The single nanocluster Hamil- 
tonians are 



% = -^ B • S, - D{Stf + D 2 (Sf) 2 + D A [(S?) 2 (S?) 2 + (5?) 2 (5?) 2 + (S?) 2 (Sf) 2 ] , (17) 
with the total magnetic field 

B = B e z + He x . (18) 
The interaction term takes into account the dipolar spin interactions 

Ha = Y. D t S ? S ^ (19) 

a/3 

through the dipolar tensor D°f = jj^ {d a p — 3ngn^j /r 3 -, in which = | r^- 1 , = r^/r^, 
and r y = r j — r j . 

One sometimes says that spin systems are similar to atomic systems, where transition 
dipoles are correlated by means of the photon exchange through the common radiation field. 
This correlation leads to coherent atomic radiation called the Dicke superradiance [34] . One 
says that moving spins also radiate electromagnetic field that could yield the correlated spin 
motion, in the same way as in the Dicke effect. To check whether this is so, we need to 
compare the time t ra d, required for inducing spin correlations through the common radiation 
field with the spin dephasing time T 2 . As the radiation time [35,36] for nanoclusters, we 
have 

3c 3 

trad — ofc 2 SO ' V^O) 

while the spin dephasing time is 

For the typical nanocluster density p ~ 10 20 cm" 3 and S ~ 10 3 , the spin dephasing time is 
T 2 ~ 1(T 10 s. While for the radiation time (20), with uj ~ 10 11 s _1 , we have t rad ~ 10 8 s = 10 
years. The ratio of times (20) and (21) is extremely large: t ra d/T 2 = 3c 3 p/(2u 3 ) ~ 10 18 . This 
tells us that the spin motion in no way can be correlated through electromagnetic radiation. 
That is, the Dicke effect has no relation to the coherent spin motion. But spins can be 
correlated only through the Purcell effect requiring the presence of a feedback field caused 
by a resonator. 

The feedback rate due to the resonator is 

7o = TiVP^llS . (22) 

The reversal time for iV correlated nanoclusters becomes 

1 T/ + 1 

. — __ — v coil _ b rev (O^ 

rev ~ 70 " ttH^SN - N ' [ 6) 

where t\ ev is the relaxation time (15) for a single nanocluster inside the same coil. 

We solved the evolution equations for the nanocluster assemblies involving the scale 
separation approach [4,14] that is a generalization of the Krylov-Bogolubov [37] averaging 
method. Four classes of spin objects have been investigated. 
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(i) Polarized nuclear materials, such as propanediol C3H8O2, butanol C4H9OH, and am- 
monia NH 3 , with the parameters: S = 1/2, p = 10 22 cm" 3 ,T = 0.1K,S ~ 10 4 G,Wo ~ 
10 V 1 , A ~ 10 2 cm,T! ~ 10 5 s,T 2 ~ 10~ 5 s, r = I/7 ~ 10" 6 s . Recall that in these 
nanomolecules the magnetization is due to polarized proton spins. 

The following characteristic times are found: thermal-noise time t T ~ 10 16 s ~ 10 9 years, 
radiation time t ra( i ~ 10 15 s ~ 10 8 years, and reversal time t rev ~ 10 _6 s. 

Therefore, neither the Nyquist thermal noise nor the photon exchange through the radi- 
ated field play any role in the relaxation process. Spin dynamics, resulting in the magnetiza- 
tion reversal, is completely due to the action of the resonator feedback field. As is explained 
above, the same concerns nanomolecules and nanoclusters 

(ii) Nuclear polarized ferromagnets, where proton spins are polarized and interact through 
hyperfine forces with electrons participating in forming ferromagnetic order. In such ma- 
terials, the electron subsystem plays the role of an additional resonator enhancing effective 
nuclear correlations. Being interested in the motion of nuclear spins, under a fixed mean 
electron magnetization, we find the reversal time t rev ~ 1CT 9 s. 

(iii) Molecular magnets, such as Mni2 and Fe§, with the typical parameters: S = 10, p = 
10 2o _ io 21 cm - 3 ,T B = IK, Bq ~ 10 5 G,w ~ 10 13 s _1 ,A ~ 10- 2 cm,w A = E A /h ~ 10 10 - 
10 12 s -1 ,7i ~ 10 5 - 10 7 s,T 2 ~ 10~ 10 s. The reversal time is t rev ~ 10~ n s. 

(iv) Magnetic nanoclusters composed of Fe, Ni, and Co, at T = 1 K, with the typical 
parameters: S = 10 3 ,p = 10 20 cm- 3 ,T B = 10 - 40K, B ~ 10 4 G = lT,w ~ lO 1 ^ -1 ^ ~ 
lcm, Ti ~ 10 34 s ~ 10 27 years,T 2 ~ 10~ 10 s. The reversal time can be very small reaching the 
value t rev ~ 10~ 12 s. 

In the case of magnetic molecules and, especially, nanoclusters, because of their high 
spins, the system of many clusters can produce quite strong coherent radiation of the maximal 
intensity 

I max - |f S^Nl h , (24) 

where N co h ~ pA 3 is the number of clusters in a coherent packet. The intensity of radiation 
of magnetic molecules, with N co h ~ 10 14 , is of order I max ~ 10 5 Vr. And for magnetic 
nanoclusters, with N co h ~ 10 20 , the radiation intensity can reach I max ~ 10 12 VF. 

There can happen several regimes of spin dynamics depending on the initial spin polar- 
ization, the strength of a triggering pulse, and the effective coupling parameter 

where 72 = 1/T 2 and A = 00 — Uq is the detuning from resonance. These regimes for 
nanoclusters can be classified analogously to those occurring for nuclear magnets [4,16]: 
incoherent free relaxation, weakly coherent free induction, weakly coherent superradiance, 
strongly coherent pure superradiance, strongly coherent triggered superradiance, pulsing 
superradiance, and punctuated superradiance [38]. 

It is important to stress that the existence of magnetic anisotropy in magnetic nanoclus- 
ters does not preclude the realization of fast spin reversal, provided the external magnetic 
field is sufficiently strong. The influence of the anisotropy energy Ea on the spin reversal 
of a nanocluster system is shown in Fig. 2, where A = Ea/ojq. This regime corresponds to 
pure spin superradiance. 
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Figure 2: Influence of magnetic anisotropy on the 
spin reversal of an ensemble of many nanoclusters, 
with the parameters typical of Fe, Ni, and Co. 
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Coherent dynamics in the spin assemblies, formed by magnetic nanomolecules, have an 
important difference from the spin dynamics in an ensemble of magnetic nanoclusters. Mag- 
netic molecules are identical and form the systems with well organized crystalline lattices. 
While magnetic nanoclusters vary in their shapes, sizes, and total spins, which results in an 
essential nonuniform broadening. Computer simulations, accomplished together with V.K. 
Henner and P.V. Kharebov, demonstrate that this nonuniformity does not destroy coherent 
spin motion. A detailed analysis of the computer simulations, with nonuniform nanocluster 
distributions, will be presented in a separate publication. 

In conclusion, we have considered spin dynamics in magnetic nanomolecules and nan- 
oclusters, starting from a strongly nonequilibrium state, with the magnetization directed 
opposite to the applied external magnetic field. We have compared several methods of de- 
scribing the spin dynamics, showing that a microscopic approach, based on the quantum 
equations of motion, is the most accurate. We also have analyzed the influence of different 
effects on spin dynamics. The effects of the Nyquist-noise triggering and of Dicke correla- 
tions are found to be negligible for spin systems. This principally distinguishes spin systems 
from atomic systems or quantum dot systems [39], where correlations, leading to coherent 
radiance, are caused by the Dicke effect of interactions through the common radiation field. 

The feedback field, developing in the resonator, reaches rather high values, of the order 
of the applied constant magnetic field. Such a strong feedback field suppresses the influence 
of mutual cluster interactions. Generally, in an ensemble of nanoclusters of sufficiently high 
density, in addition to dipole interactions, there can appear exchange interactions [40] that 
can influence equilibrium properties of nanoclusters. But in the considered case of strongly 
nonequilibrium spin dynamics, the exchange interactions are also suppressed by the self- 
organized resonator feedback field. 

This important conclusion can be formulated as follows: Coherent spin dynamics are 
completely governed by the Purcell effect that is caused by the action of the resonator feedback 
field. 

References 

[1] Barbara B, Thomas L, Lionti F, Chiorescu I, and Sulpice A 1999 J. Magn. Magn. Mater. 
200 167 

[2] Wernsdorfer W 2001 Adv. Chem. Phys. 118 99 



9 



[3] Ferre J 2002 Topics Appl. Phys. 83 127 

[4] Yukalov V I and Yukalova E P 2004 Phys. Part. Nucl. 35 348 

[5] Wang J and Zeng X C 2009 in: Nanoscale Magnetic Materials and Applications 
(Springer, Berlin), p. 35 

[6] Bedanta S and Kleemann W 2009 J. Phys. D 42 013001 

[7] Berry C C 2009 J. Phys. D 42 224003 

[8] Beveridge J S, Stephens J R, and Williams M E 2011 Annu. Rev. Anal. Chem. 4 251 

[9] Hoang V V and Ganguli D 2012 Phys. Rep. 518 81 

[10] Bauer M, Fassbender J, Hillebrands B, and Stamp R L 2000 Phys. Rev. B 61 3410 

[11] Yukalov V I and Yukalova E P 2012 J. Appl. Phys. Ill 023911 

[12] Purcell E M 1946 Phys. Rev. 69 681 

[13] Bloembergen N and Pound R V 1954 Phys. Rev. 95 8 

[14] Yukalov V I and Yukalova E P 2000 Phys. Part. Nucl. 31 561 

[15] Yukalov V I 2002 Laser Phys. 12 1089 

[16] Yukalov V I 2005 Phys. Rev. B 71 184432 

[17] Yukalov V I, Henner V K, Kharebov P V and Yukalova E P 2008 Laser Phys. Lett. 5 
887 

[18] Yukalov V I, Henner V K and Kharebov P V 2008 Phys. Rev. B 77 134427 

[19] Yukalov V I and Yukalova E P 2011 Laser Phys. Lett. 8 804 

[20] Landau L D and Lifshitz E M 1935 Phys. Zeits. Sow. 8 153 

[21] Henner V, Raikher Y and Kharebov P 2011 Phys. Rev. B 84 144412 

[22] Gilbert T L 1956 PhD Thesis (Illinois Institute of Technology, Chicago) 

[23] Bloch F 1946 Phys Rev. 70 460 

[24] Abragam A and Goldman M 1982 Nuclear Magnetism: Order and Disorder (Clarendon, 
Oxford) 

[25] Yukalov V I 1992 Laser Phys. 2 559 

[26] Fokina N P, Khutsishvili K O and Chkhaidze S G 1992 J. Exp. Theor. Phys. 102 1013 

[27] Khutsishvili K O and Chkhaidze S G 1992 Physica B 176 54 

[28] Fokina N P, Khutsishvili K O and Chkhaidze S G 1992 Physica B 179 171 

10 



[29] Yukalov V I 1996 Nucl. Instrum. Methods Phys. Res. A 37 345 

[30] Yoshimi A, Asahi K, Sakai K, Tsuda M, Yogo K, Ogawa H, Suzuki T and Nakagura M 
2002 Phys. Lett. A 304 13 

[31] Yukalov V I 1995 Phys. Rev. Lett. 75 3000 

[32] Yukalov V I 1996 Phys. Rev. B 53 9232 

[33] Rueckriegel A, Kreisel A and Kopietz P 2012 Phys. Rev. B 85 054422 
[34] Dicke R H 1954 Phys. Rev. 93 99 

[35] Yukalov V I and Yukalova E P 2005 Laser Phys. Lett. 2 302 
[36] Yukalov V I 2005 Laser Phys. Lett. 2 356 

[37] Bogolubov N N and Mitropolsky Y A 1961 Asymptotic Methods in the Theory of Non- 
linear Oscillations (Gordon and Breach, New York) 

[38] Yukalov V I and Yukalova E P 2002 Phys. Rev. Lett. 88 257601 

[39] Yukalov V I and Yukalova E P 2010 Phys. Rev. B 81 075308 

[40] Kondratyev V N and Lutz H O 1998 Phys. Rev. Lett. 81 4508 



11 



